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ABSTRACT 

The rise in the cosmic-ray positron fraction measured by the PAMELA satellite is 
likely due to the presence of astrophysical sources of positrons, e.g. pulsars, on the 
kpc scale around the Earth. Nevertheless, assessing the properties of these sources 
from the positron data requires a good knowledge of the secondary positron com- 
ponent generated by the interaction of cosmic rays with the interstellar gas. In this 
paper, we investigate the impact of the spectral hardening in the cosmic-ray proton 
and helium fluxes recently reported by the ATIC2 and CREAM balloon experiments, 
on the predictions of the secondary positron flux. We show that the effect is not neg- 
ligible, leading to an increase of the secondary positron flux by up to ~60% above 
~100 GeV. We provide fitting formulae that allow a straightforward utilization of 
our results, which can help in deriving constraints on one's favorite primary positron 
source, e.g. pulsars or dark matter. 
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1 INTRODUCTION 



The increase in the positron fraction abov e a few GeV 
repor ted by the PAMELA collaboration |Adriani et al.l 
120091 ) has triggered a lot of interpretation attempts, 
but it is now likely that it is due to positrons origi- 
nating fro m conventiona l astrophysic a l sources, like pul- 



sars (e.g. ISturro"ck1 Il970l; jh e"nlll970l; lHarding fc Ramatyl 



1987i;lBoulares|ll989l;lAharonian et aT1ll995l; IChi et all 
Profumol 12008: iMalvshev et alj 12009; lYiiksel et al 
Delahave et alj 1201(f ) or SNRs (e.g. iBerezhko et al 



199C 



200£ 



2003 



Blasill2009l ). while som e other spatia l or so l ar effects migh t 



also play a role (e.g. IShaviv et all 120091 : iRobertsI 12010 ) . 

These works have notably shown that a very few sources 
may dominate the high energy positron flux at the Earth, 
opening interesting perspect ives for more accura te predic- 
tions in the near future (see iDelahave et al. 2010., for a de- 
tailed analysis). Nevertheless, these perspectives are based 
on the assumption that the secondary positron flux predic- 
tion is under control. 

Secondary positrons originate from nuclear interac- 
tions of cosmic-ray nuclei with the interstellar gas, and 
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have b een investigated into det ails by Moskalenko fc Strond 
(1998), and more recently by IDelahave et al.l (|2009l ). who 
provided more insig hts on the theoretical uncertainties. 
IDelahave et al.l (|2010l ) have improved the predictions of the 
latter by including the Klein-Nishina corrections to the en- 
ergy loss treatment. All these predictions rely on cosmic- 
ray nuclei spectra constrained from rather low energy data 
(< 100 GeV), which are mostly p ower laws. Neverth eless, 
two ballo on experiments, ATIC2 (|Panov et al.1 120091 ) and 
CREAM (|Ahn et al.ll2010l ) , have recently reported on a clear 
and almost concordant hardening in these spectra around 
a few TeV, with a very good statistics. Since high en- 
ergy stable nuclei have quite long range propagation (e.g. 
iTaillet fc~M aurin 2003), it is likely that this spectral inflec- 
tion is not merely local and pertains over a few kpc scale 
around the Earth. Therefore, a consistent prediction of the 
secondary positron flux should take it into account. 



In this paper, we study into detail the impact of these 
new cosmic ray measurements on the secondary positron flux 
predictions. We provide the reader with user-friendly fitting 
formulae which summarize our results, and which can be 
used to constrain any extra source of cosmic-ray positrons. 
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2 FROM THE HARDENING OF COSMIC-RAY 
NUCLEI SPECTRA TO THE HARDENING 
OF THE SECONDARY POSITRON 
SPECTRUM 

2.1 Generalities 

Predictions of secondary positrons are usually valid in 
the frame of a propagation model, which specifies the 
way the propagation equation is solved. For more in- 
sights o n cosmic ray propagat i on, w e refer the reader 
to e.g . iGinzburg fc Svrovatskiil (|l964h : | Berezinskii et all 
(1990); Longai r) dl994Tl. Here, w e adopt the formalism de- 
scribed in lDelahave et aD ()2010l ). where convection and dif- 
fusive reacceleration are neglected, which is known to be a 
pod approximation i n the GeV-TeV energy energy range 
Delahave et al.ll2008l) . and where fully relativistic energy 
losses are considered. Aside from energy losses, our propaga- 
tion ingredients are therefore the diffusion coefficient K(E), 
that we take homogeneous, and the half-thickness L of the 
diffusion slab. We use sets of propagation parameters con- 
sistent with the analysis of the secon d ary-to -primary nu- 
clei rat ios performed bv IMaurin et al.l (|200lh . and further 
used in iDonato et alj (|2004l ) to define the minimal, median 
and maximal sets widely used in the literature. We note 
th at the med model, w hich was the best-fit model found 
in IMaurin et al. 

I (|200ll ). has properties very similar to the 
best-fit model derived more recently in IPutze et al.l (|2010h 
from a Markov chain Monte Carlo analysis. Details on the ef- 
fects of these param eters on the electron an d positron prop- 
agation are given in lDelahave et all (|2010l ). For the energy 
losses, we adopt the model denoted Ml in this reference. 

Our whole propagation framework can be encoded in 
the form of a Green function Q(E,xq <— E s ,x s ) that char- 
acterizes the probability of a positron injected at any coor- 
dinate x s with energy E s to reach an observer on Earth with 
energy E ^ E s . This allows us to write the positron flux as 
the following convolution: 



4>{E) 



d^Xt 



dE s g(E,x Q <- E s ,Xs)Q{E 3 ,x s 



(1) 



where Q is the source term that we are going to deter- 
mine in the following . For secondaries, it formally reads (e.g. 
iDelahave et~aii r2009) 



Q(E,x) = 4tt^ J dE k <j>i{E k ) ^§-(E k -> E)n 3 (x) 



(2) 



where (j>i is the flux of a cosmic-ray species of index i, rij 
is the interstellar density of a gas species of index j, and 
doij is the inclusive nuclear cross section associated with 
the production of a positron of energy E from an ion of 
kinetic energy E k . For the cosmic-ray nuclei and interstel- 
lar gas, we can safely consider the dominant species only, 
i.e. the protons and a ions on the one hand, and the hy- 
droge n (90%) and h e lium (10%) gas on the other hand. 
As in IDelahave et~ai1 l|2009l ) we assume an overall gas den- 
sity of no = 1 cm -3 homogeneously distributed inside an 
infinite flat disk of half-thickness h = 100 pc, such that 
ntot(x) = 2 hno S(z), where z is the coordinate perpendic- 
ular to the Galactic plane. These values are justified either 
by measurements of the interstellar medium (|Ferrierej|200ll ) 
and the fact that high energy positrons have very short range 
propagation due to efficient energy losses — large scale fluc- 



tuations of the gas density have almost no effect on the local 
high energy positron flux. 



2.2 The incident cosmic ray flux 

In IDelahave et all (Hp09), we solved Eqs. CD and © as- 
sumin g the proton and a fluxes as fitt ed in IShikaze et al.l 
(120071) to the low energy BESS da ta (|Sanuki et al.l |2000| ; 
I Wang et al1l2002l ; lHaino et al.ll2004r ). These fits are recalled 
hereafter: 



s (E k ) = Ap ai 



n 



<t>a SS (E k /n) 



BP" 



1GV 



1GV 



(3) 



= (1.94 cm -2 s~ 
= (0.71 cm -2 s _1 
laws of indices 



GeV 
GeV" 
a2 and 



1 sr" 1 , 0.7, 2.76) 
sr" 1 , 0.5, 2.78), 
&2, respectively, 



with (A,ai,a2) 
and (f?,bi,&2) : 
i.e. single power 
at high energy. 

These functions are displaye d (solid curves 
Fig. [D t o gether with the BE S S JSanuki et al 



Wan g et all 2002: 



Wang 
2007fl . 



CAPRICE 



lHaino et al] 12004 
dBoezio et al. 



rrves) 
2000: 



.Shikaze et al.l 
19991. l200liT 



AMS 01 (|Aguilar et aD |2002|). ATIC2 l|Panov et all I2006L 
120091) , and CREAM (|Ahnetai1 \20ld ) data — note that 



these data have been corrected for solar modulation effects 
(demodu lated) by means of th e Force Fi eld approxi- 
mation ijGleeson fc Axfordl 1 19681 : iFiskl Il97ll ). with Fisk 
potentials made explicit in the plot. It is clear that though 
these parameterizations provide reasonably good fits to the 
low energy data, they completely fail above a few tens of 
GeV. In particular, we can remark that the proton data 
(left panel) are overshot between ~ 10 GeV and a few TeV, 
while the helium data (right panel) are underpredicted 
above ~ 100 GeV/n. Moreover, though CREAM and 
ATIC2 seem to agree in their measurements of the helium 
flux, there is an unequivocal discrepancy in the proton flux. 

Since the f unct ions of Eq. (13) were used 
in IDelahave et all (120091 ) and IDelahave et all (|2010T ) to 
improve the predictions of the secondary positron flux, 
it is worth revisiting these predictions again in light of 
the new cosmic ray data. For the inclusive nuclear cross 
sec tions, we stil l use the numerical approach presented 
in iKamae et al.l (|2006T ) for the pr oton-proton collision, 
with the correction prescriptions of iNorburv fc Townsendl 
|2007h for nucleus-nucleus interactions. The change in the 
derivation of the positron source term defined in Eq. ([2]) 
will therefore only come from the updated fits of the proton 
and helium fluxes. 

To illustrate the difference in the positron source term 
arising from considering either the CREAM or the ATIC2 
proton data, we consider two different modelings defined 
by the preference we put in one or the other experiment. 
The proton flux parameterization will be the only change 
between these two cases, since the helium data of the two 
experiments are consistent and can be simultaneously fitted 
by the same function. 

In the following, we denote Flp and F2p the parameter- 
izations associated with the proton data, the former provid- 
ing a good fit to the CREAM data, and the latter providing 
a good fit to the ATIC2 data. We also define a function FIHe 
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Figure 1. Left panel: cosmic-ray proton data. Right panel: cosmic-ray a data. 



that provides a good fit to both CREAM and ATIC2 helium 
data. Functions Flp, F2p and FIHe are given in Sec. 13.11 
and are displayed in Fig. [TJ in the left panel for the two 
first, and in the right panel for the last one. The functions 
of Eq. (|3]) will be referred to as low energy fit, or reference 
fit. In contrast, the combination of Flp and FIHe will be 
referred to as CREAM fit, while F2p and FHe as ATIC2 fit. 



2.3 Impact of the cosmic ray hardening on the 
secondary positron source term 

In the left panel of Fig. [2j we derive the secondary positron 
source term defined in Eq. ([2]) and associated with the low 
energy fit (black curve), the CREAM fit (red curve), and the 
ATIC2 fit (blue curve); note that all of them run close to 
a power law in energy of index -2.7. The contribution com- 
ing from the cosmic-ray proton is shown explicitly (dashed 
curves) for each modeling, as well as the a contribution (dot- 
ted curves). We see that the relative increase in the a con- 
tribution, from the low energy fit to FIHe, is very large, 
reaching a factor of ~ 3 around a few TeV. Nevertheless, 
it remains subdominant with respect to the proton contri- 
bution. The difference in the proton flux modeling trans- 
lates almost linearly into the positron source term, lead- 
ing to a relative decrease with respect to the low energy 
fit below 100 GeV (1 TeV), and a relative increase above, 
for the ATIC2 (CREAM, respectively) fit. This comes from 
the fact that the in clusive cross section f eaturing in Eq. (J5J) 
scales like ~ 1/Ek (iDelahave et al.ll2009f) . straightforwardly 
leading to Q tx cj>i(Ek); this explains why all contributions 
almost scale like E~ 2 ' 7 . When summing up the contribu- 
tions coming from proton and a interactions, we see that 
the relative decrease apparent in the proton-only case is less 
prominent due to the positive yield, though modest, of the a 
interactions: the net effect is a very slight decrease below 100 
GeV, and a larger increase above. This is illustrated in more 
detail in the right panel of Fig. [2] where we plot the relative 



difference of the CREAM (dashed curve) and ATIC2 (dot- 
ted curve) cosmic ray-induced positron source terms with 
the low energy reference case: the slight decrease below 100 
GeV makes a 10% difference at most with the reference case, 
while the increase above reaches ~ 30% (60%) above a few 
TeV for the CREAM (ATIC2, respectively) configuration. 
The impact of using these new cosmic ray data is therefore 
not negligible in terms of secondary positron production. 



2.4 Updated predictions for the secondary 
positron flux 

Because of propagation, guessing the effect of the new cos- 
mic ray modelings on the secondary positron flux prediction 
might not look, a priori, as straightforward as guessing their 
effect on the positron source term. This is formally due to 
the non-trivial dependence of the positron Green fun ction on 
energy. We refer the reader to lDelahave et all (|201Ch . in par- 
ticular to their Sect. 2.3, for more insights about this depen- 
dence. Nevertheless, since here we deal with a source term 
which is homogeneously distributed inside a thin disk of half- 
thickness h, and with an energy dependence close to a power 
law, we can assume that Q(E,x m 2 h Qo S(z) (E/GeV)' 1 . 
In that ca se, the positron flux a t the Earth can be approxi- 
mated as (jDelahave et al.ll2010T ): 

(£)£^L%(£/GeVH, (4) 

where Ko is the normalization of the diffusion coefficient, t; 
is the energy loss timescale, and 7 « 7 + 0.5 (a + 8 — 1) is the 
predicted flux spectral index which depends on the injected 
index 7 and also on the diffusion coefficient index 8 and the 
energy loss index a (a — 2 in the Thomson approximation, 
but is < 2 when Klein-Nishina corrections become sizable). 

The above equation is of great interest to anticipate the 
coming results, since it tells us that the ratio of two differ- 
ent flux predictions scales like the ratio of the corresponding 
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Figure 2. Left panel: Secondary positron source term for three different fit-based assumptions for the incident cosmic ray spectra; the 
contributions due to cosmic-ray protons (dashed curves) and a particles (dotted curves) are shown explicitly aside from the overall 
contributions (solid curves). Right panel: relative difference between the obtained overall source terms. 



source terms at any energy. Therefore, the relative differ- 
ences in flux predictions should be very close to the relative 
differences in the corresponding source terms, the latter be- 
ing already plotted in the right panel of Fig. [2] This is ac- 
tually verified in the right panel of Fig. O which can hardly 
be differentiated from the one previously mentioned, illus- 
trating how efficient the approximation given in Eq. Q is in 
this context. In the left panel of Fig. [3] we show the corre- 
sponding predictions for the secondary positron flux, which 
are derived with the med propagation parameters and with 
model Ml for the energy losse s — these parameters can be 
found in lDelahave et alT(|2010r ). From these plots, we readily 
conclude that the secondary positron flux prediction is af- 
fected over the whole GeV-TeV energy range, decreasing by 
~ 10% below ~ 100 GeV, and increasing by more than 30% 
at TeV energies, up to 60% in the case of the ATIC2 cosmic 
ray fit. The spectral index is therefore increased accordingly 
by almost 2 digits, from ~-3.5 to ~-3.3 in the med case em- 
ployed here. 

For the sake of completeness, it is interesting to 
derive the theoretical uncertainty bands associated with 
these novel predictions, which come from the uncertain- 
ties in the p ropagation parameters . To this aim, we 
proceed as in IDelahave etafl (|2009ft by bracketing the 
med with the min and max propagation configurations. 
We report our results for the positron flux (fraction) 
in the top (bottom, respectively) panels of Fig. [4] — a 
solar modulat ion is applied, us i ng the For c e Fie ld ap- 
proximation l|Gleeson fc Axfordl 1 19681 ; iFiskl llQTlh with 
a Fisk potential of 600 MV, and ignoring any charge 



the positron fraction / = 4> e +/(4' e + + <f> e -)i we deter- 
mine d the denominator above 10 GeV from the Fermi -LAT 
data (|Abdo et al . 20091; iFermi LAT collaboratiorj|20icft . and 



GeV (e.g. iHeber fe Potgiete: 



ly rnip c 
rl 12006ft 



T he positron flux 



data a r e taken from Boezio et al. ( 2000ft ; iDuVernois et al] 
(2001ft ; lAguilar et all (120021) while the positron f ractio n 
data are from Barwick et alj (Il997rj ; Beattv et all l|2004ft ; 
lAguilar et "all (|2007r i; lAdriani et alj l|2010ft . Tocalculate 



used the AMS01 data (Aguilar et alj 20021 ) to constrain the 
electron flux at lower energy. In all panels, we also display 
the re sult obtained by iMoskalenko fc Strong] (|l998T ). as fit- 
ted bv lBaltz fc Edsiol j 1998ft . since it is still very often used 
in the literature as a reference. Note that the predictions 
obtained from the low energy cosmic ray fit a r e obv iously 
identical to the ones derived in IDelahave etafl (|2010ft . 

From the predictions based on the low energy cosmic ray 
fit to the ones based on the ATIC2 fit, the secondary positron 
spectrum is hardened, which translates into a slightly flatter 
positron fraction. Of course, we did not expect the present 
study to be relevant to the discussion on the rise of the 
positron fraction itself, since the enhancement in the sec- 
ondary positron flux was already known to be by far too 
small from the incident cosmic-ray nuclei data. It is instead 
very useful so as to constrain any extra source of positrons, 
like pulsars or dark matter, with the data. Indeed, in that 
case, one needs to add a secondary contribution to the pri- 
mary one in a consistent manner before comparing the sum 
to the data. From the top panels of Fig. 2] we note inci- 
dentally that the min propagation setup already leads to a 
conflict with the data because of its too small value of Ko 
[see Eq. @] (associated with a small value of L, Ko/L be- 
ing roughly fixed by secondary-to-primary cosmic-ray nuclei 
ratios). This configuration is in any case obsolete, since it 
is no long er supported by re cent secondary-to-primary anal- 
yses (e.g. iPutze et alj|20icft. nor by the Fermi-LAT diffuse 
gamma-ray data ( e.g. lAbdo et alj|2010ft (see also a dedi- 
cated discussion in iLavalld 120101 ) . Nevertheless, it can still 
be thought of as an extreme configuration and used for il- 
lustration purposes. 

We provide user-friendly empirical fitting functions as- 
sociated with all these propagation models in Sec. 13.21 
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Figure 3. Left panel: Secondary positron flux predictions associated with the 3 different fit-based assumptions considered for the incident 
cosmic ray spectra; the mecipropagation setup has been used. Right panel: relative difference between the predictions shown in the left 
panel. 





Figure 4. Top panels: secondary positron flux predictions, both interstellar (IS) and on top of atmosphere (TOA), assuming a Fisk 
potential of 600 MV for the solar modulation. Bottom panels: associated positron fractions. These predictions correspond to source terms 
calculated from cither the low energy cosmic ray fit (left panels), the CREAM cosmic ray fit (middle panels), or the ATIC2 cosmic ray 
fit (right panels). 



6 J. Lavalle 



3 USER-FRIENDLY FITTING FORMULAE 



4 CONCLUSION 



In this section, we provide the reader with the parametric 
functions (i) that we used to fit the cosmic ray data and 
(ii) also the ones which fit our prediction results for the 
secondary positron flux. 



3.1 Cosmic-ray proton and helium interstellar 
fluxes 

The proton data can be accommodated by two different 
empirical functions, depending on whether one favors the 
CREAM (Flp) or ATIC2 (F2p) data at high energy. These 
functions are both characterized by the following parame- 
terization: 



Fp(£ fc ) = 



E k 

Epl, 



lOGeV/n 



-~ip 



1 + 



E p2 



1 + 



Ek_ 

E pa 



(5) 



One recognizes a standard power law of main index y p asso- 
ciated with an exponential attenuation factor active at low 
energy, which clearly improves the low energy fit with re- 
spect to Eq. (J5J, and a double spectral correction of indices 
P2 and P3 operating above kinetic energies E P 2 and E p s, re- 
spectively. Note that function Fp behaves asymptotically as 
a power law of index = —7 + p 2 + Pa- 
in contrast, a single function (FIHe) is enough in the 
case of helium ions, since both sets of high energy data are 
in agreement. In that case, a slight correction to Eq. ([3]) is 
enough, so that 



FHe(£ fc /n) = <f> b a eBB (E k /n) 



TZ 

K p2 



1 + 



TZ 

TZ p3 



(6) 



where <^>a ess is given in Eq. ((3]), and TZ is the rigidity. There 
is also a double spectral correction as in the case of protons. 
As function Fp, function FHe behaves asymptotically as a 
power law of index 700 = — 62 + P2 + P3, where 62 is defined 
in Eq. ©. 

Given the spread in the available data, we did not per- 
formed a x 2 selection of the parameters. Indeed, this would 
require an a priori or expert selection of the data, which is 
beyond the scope of this paper. The values of the parameters 
used for functions Flp, F2p and FIHe are listed in Table [T] 



3.2 Secondary positron flux predictions 

Here, we define an empirical function that provides a very 
good fit to the interstellar secondary positron flux predic- 
tions presented in Sec. 12.41 



'fit (E) = exp <{ ^2 a 

. i=0 



In 



— ) 

GeV J 



(7) 



The parameters associated with all the configurations dis- 
cussed in Sec. 12.41 are given in Table f2] We emphasize that 
these parameters are valid only for the propagation models, 
energy loss configurations, and incident cosmic ray model- 
ings discussed through this paper. These empirical fitting 
functions should not be used for other sets of parameters. 
They are valid from ~ 0.1 GeV to ~ 10 TeV. 



In this paper, we have studied the impact of the spectral 
hardening observed in the cosmic-ray proto n and helium 
fluxes by the ATI C2 (|Panov et alj|2006l . l200gh and CREAM 
I Ahn et alj I2010T ) experiments, on the secondary positron 
flux prediction. To this aim, we have revisited the calcula- 
tion of the secondary positron source term, which spatially 
tracks the interstellar gas, showing that its energy distribu- 
tion roughly scales like the incident cosmic ray spectrum. 
Because of the discrepancy in the proton fluxes observed 
by the CREAM and ATIC2 balloons — in contrast, both 
agree on the helium flux — we have considered two different 
cosmic ray modelings: one based on a fit to the CREAM 
proton data (moderate case), and another based on a fit to 
the ATIC2 proton data (maximal case), both using the same 
helium flux parameterization. The former (latter) case led 
to a 30% (60%, respectively) increase in the production rate 
of TeV positrons. 

Then, we have propagated these positrons to 



the Earth, 



using 



in iDelahave et al 



the 



(201 



0). 



ropagation framework described 
which includes a fully relativis- 
tic treatment of the energ y losses, with spatial diffusion pa- 
rameters as constrained in lMaurin et al.l (|200ll ) , still consis- 
tent w ith the more recent analysis performed in lPutze et al.l 
(2010). We have notably explained why, in this context, 
the relative differences in the differential flux predictions 
were almost equal to the relative differences in the energy- 
dependent positron production rate, and consequently to 
the relative differences in the considered incident cosmic ray 
fluxes. This led us to establish a robust estimate of the effect: 
as for the positron production rate, the secondary positron 
flux is increased by up to 30% (60%) at TeV energies when 
constraining the cosmic-ray nuclei fluxes from the CREAM 
(ATIC2, respectively) data. Theref ore, these predictions 



differs from the ones p erform ed in Moskalcnk o Sz Strond 
(1998), IDelahave etafl ((2009) and IDelahave et al. 

I (|2010h . 

resulting in harder secondary positron spectra. We have 
also derived an estimate of the theoretical uncertainties. 
T hese results are comple mentary to the recent analysis 
of lDonato fc Serpicol (|2010h on the secondary antiproton and 
diffuse gamma-ray fluxes. 

Finally, we have provided the reader with user-friendly 
parametric functions that allow to reproduce all the results 
derived in this paper, so that they can be straightforwardly 
exploited. For instance, one can use these functions to con- 
strain any extra source of primary positrons against the ex- 
isting (or forthcoming) data in a self-consistent way — in 
which case one needs to sum up the secondary and primary 
contributions. We again emphasize that these results are 
valid only in the frame of the propagation models described 
through the paper, so if used, one has to make sure that the 
same parameters are employed for the primary component. 
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[TeV] 
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Flp 


3.09 x 1CT 3 


2.8 


4 
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-0.29 
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2.8 
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1.05 


1.5 


0.4 


10 


-0.35 
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1 


0.5 


10 


-0.5 



Table 1. Parameters used to fit the cosmic-ray proton (Flp, F2p) and helium (FIHe) data, according to Eqs. l(5} and JSJl. For function 
FIHe, any parameter E p (in GeV) in this table corresponds to 1Z P (in GV) in the associated equation. 







Low E CRs 






CREAM CRs 






ATIC2 CRs 






min 
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med 


max 
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-4.61 
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-6.4 


-4.61 
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-6.4 


-4.41 


-5.48 


-6.4 
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-3.55 


-3.486 


-3.37 


-3.6 


-3.53 


-3.41 


-3.6 


-3.53 


-3.41 


c 2 (xl0" 2 ) 
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-8.2 


-9.28 


-8.99 


-8.83 


-8.98 


-8.68 


-8.52 


c 3 (xl0- 2 ) 
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2.13 


2.73 


2.67 


2.64 


2.79 


2.73 


2.69 


c 4 (xlO- 3 ) 


-1.41 


-1.38 


-1.37 


-1.8 


-1.77 


-1.754 


-1.88 


-1.84 


-1.82 



Table 2. Parameters used to fit the interstellar secondary positron flux predictions, from Eq. (0. 
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